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Abstract 



We study a c = — 2 conformal field theory coupled to two-dimensional quantum 
gravity by means of dynamical triangulations. We define the geodesic distance r 
on the triangulated surface with N triangles, and show that dim[r H ] = dim [AT], 
where the fractal dimension du = 3.58 ± 0.04. This result lends support to the 
conjecture du = — 2oti/ot_i, where a_ n is the gravitational dressing exponent of 
a spin-less primary field of conformal weight (n + 1, n + 1), and it disfavors the 
alternative prediction du = — 2/y s tr- On the other hand, we find dim[/] = dim[r 2 ] 
with good accuracy, where I is the length of one of the boundaries of a circle with 
(geodesic) radius r, i.e. the length I has an anomalous dimension relative to the area 
of the surface. It is further shown that the spectral dimension d s = 1.980 ±0.014 for 
the ensemble of (triangulated) manifolds used. The results are derived using finite 
size scaling and a very efficient recursive sampling technique known previously to 
work well for c = —2. 
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1 Introduction 



Liouville theory and matrix models have been successful in explaining a number 
of features of conformal field theories coupled to two-dimensional quantum gravity. 
However, our primary interest in a theory of quantum gravity concerns geometry: 
which concepts of geometry survive the quantum average, and how is the geometry 
changed by this average. These questions have not been clarified by Liouville field 
theory or matrix model techniques. In the last couple of years significant progress has 
been made in this direction, starting with the introduction of the so-called transfer 
matrix |1| . It was shown that a reparametrization invariant formulation of the two- 
point function in quantum gravity has a simple geometric interpretation || and that 
a generalization of the two-point function to include matter fields allows a geometric 
interpretation of the /fPZ-exponents ||, an interpretation first conjectured in [Q. 
In addition it was realized that finite size scaling analysis of the two-point functions 
were very efficient tools for extracting critical exponents |5], || . In this article we will 
take advantage of this new technology and combine it with the efficient recursive 
sampling algorithm developed earlier for c = —2 conformal field theory coupled to 
two-dimensional quantum gravity used in || where the fractal nature of quantum 
gravity in two dimensions was first numerically confirmed. 

2 The model 

The c = —2 model coupled to quantum gravity corresponds to a non-unitary (1, 2) 
conformal field theory coupled to quantum gravity, known as a topological quantum 
gravity. There exists an explicit realization of the model within the framework of 
dynamical triangulations. In this framework the partition function for c Gaussian 
fields coupled to two-dimensional quantum gravity is 

Z N = J2(detC TN y c/2 . (1) 

T N 

In ([I]) the summation is over all triangulations Tn with fixed topology (which we will 
always assume to be spherical in this paper) built from N triangles, and C? N is the so- 
called adjacency matrix of the graph corresponding to the triangulation 7/v- Notice 
that the triangulations of the spherical surface are in one-to-one correspondence 
with the 3 connected planar graphs with no external legs, and thus it is possible 
to generate any triangulation TJv if we generate and connect 3 trees and rainbow 
diagrams with the correct weight. If c = —2, the weight for generating 3 trees and 
rainbow diagrams is 1, i.e. Eq. (Q) can be written [[/j 

Z »=NT2( £ X )( £ 0' (forc=-2), (2) 

tree diagrams rainbow diagrams 

with N + 2 legs with (N + 2)/2 lines 

where the first summation is over all rooted 3 tree diagrams with N + 2 external 
legs and the second summation is over all rainbow diagrams with (N + 2)/2 lines. 
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l/(iV + 2) is the symmetry factor which comes from connecting the tree diagram 
and the rainbow diagram. Using a recursive algorithm for generating 3 trees and 
the rainbow diagrams in (||), it is possible to create a large number of independent 
triangulations T N with the weight detCV^. We refer to || and || for details. Sup- 
plementary to the first studied in |J we will in this work study the finite size aspects 
of the observables associated with Zn, the main motivation being that finite size 
scaling by far is the most reliable method for extracting continuum physics in critical 
systems and we now understand that this is true also for these systems coupled to 
quantum gravity. 

One very important point in the above setup is that we have the concept of 
distance, even if we usually associate c = —2 model with a topological gravity. To a 
triangulation TJv we can unambiguously associate a piecewise linear manifold with 
a metric dictated by the length assignment e to each link. From a practical point of 
view we use instead a graph-theoretical distance between vertices, links or triangles. 
In the limit of very large triangulations we expect that the different distances when 
used in ensemble averages will be proportional to each other. To be specific we will 
in the following operate with a "link distance" and a "triangle distance" . The link 
distance between two vertices is defined as the shortest link-path between the two 
vertices, while the triangle distance between two triangles is defined as the shortest 
path along neighboring triangles between the two triangles. In this way the triangle 
distance becomes the link distance in the dual 3 graph. 

In the following we will report on the measurement of two quantities related to 
the fractal structure of quantum space-time: The total length (I) (and the higher 
moments (/")) of boundaries of spherical balls of (geodesic) distance r, and the 
measurement of so-called spectral quantities, originating from the study of random 
walks on the manifolds. 

Let us define the observables on the triangulation 7/v more precisely. We consider 
a spherical ball of radius r and its shell for a given triangulation T^. The spherical 
ball consists of all vertices with link distance r' < r and the spherical shell consists 
of all vertices with link distance r, where the distance is measured from a given 
vertex t> which is considered as the center of the spherical ball. In the same way 
we can define the spherical shell in terms of triangle distance. We will use both 
definitions in the following. The spherical shell will in general consists of a number 
of connected components. If we take the average over all positions of vq and all 
triangulations T N , we get a distribution p N (l,r) of the length I (measured in link 
units) of the connected components of the spherical shells of radius r, i.e. 

oo 

(l n )r,N = i: in PN(l,r). (3) 
1=1 

In particular we introduce the special notation n^{r) = (l(r)) N and expect the 
fractal dimension to be related to n^{r) by 

n N {r) ~ r dH -\ 1 « r < N 1/d ' H . (4) 

According to general scaling arguments @, |, |J we expect the following behaviour 
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for njv(r): 

n N (r) ~ N 1 - 1 ^F 1 (x), x = , (5) 

and we expect Fi(x) to fall off rapidly when s>l. 

The spectral properties are derived from the study of random walks on the 
triangulated surfaces via the diffusion equation 

0M + i) = — E^K*)> ( 6 ) 

71 v (vv') 

where t is the diffusion time and the summation is over the n v neighboring vertices^ 
v' to v. We here consider the initial condition of <f)(v,t), 

(f>(v,0) = —6 VjV0 , (7) 

where v is a fixed vertex. Eqs. (H) and (|7|) are discretizations of the diffusion 
equation on a continuum manifold. We will be interested in observables associated 
with the diffusion process obtained by averaging over the chosen vertex v o as well as 
different triangulations T N . In the following discussion we will always assume that 
this average has been performed. Let us denote the probability of diffusion out to a 
link distance r in time t by k]y(r, t). By definition we have 

oo 

J2n N {r) k N {r,t) = 1 . (8) 
From scaling arguments [§, §, [L(J we expect the following scaling 

k N (r,t) = -p{x,y), x = j^, V = ^ , (9) 

where the new exponent A is defined such that y will be finite in the scaling limit. 
The spectral dimension is defined from the return probability by 

M°'*) 7^77^' for £ ~°> ( 10 ) 

while the average geodesic distance travelled by diffusion at time t is 

oo 

(r(t)) N = J2m N (r)k N (r,t)^N 1 / d "-^r, for *~0. (11) 

r=0 

If a kind of "smooth" fractal is expected here again, fcjv(0, t) and {r(t)) N exist and 
are different from zero in the limit iV — ► oo. This implies the scaling relations 



1 Also in this case we have the possibility to formulate a diffusion in terms of triangles and 
triangle distances, rather than vertices and link distances. However, the finite size effects are 
larger for triangles and we will in the following only use link distances when we discuss diffusion. 
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3 Numerical results 



3.1 The simulations 

The simulations are performed by generating a number of statistically independent 
configurations using the algorithm mentioned in the introduction (see M and pfl 



for details). We use the high quality random number generator RANLUX ||T6| , \Vq 
whose excellent statistical properties are due to its close relation to the Kolmogorov 
K-system originally proposed by Savvidy et.al.0, [I5j in 19860. We report results 



on system sizes ranging from 2000-256000 triangles. The number of configurations 
obtained depends on the lattice size and on the observable that we measure. We 
choose 20 random vertices on each configuration in order to perform correlation 
function measurements. We need to collect more statistics to test Eq. (0), where 
we have between 4.2 x 10 6 and 1.6 x 10 6 configurations. For the 128K and 256f^ 
lattices we have 6 x 10 5 and 2 x 10 5 configurations respectively. In order to mea- 
sure the moments (l n ) rN and their scaling properties we need a factor of 10 2 less 
configurations: We have approximately 50000 configurations for each lattice size. 
Unfortunately, the computer effort for making the measurements is comparable to 
the one needed to test Eq. (|5|) with enough accuracy. For the diffusion equation 
we collect 2500 configurations on which we perform 5 measurements. For the three 
largest lattices we have 620, 600 and 400 configurations respectively. 

3.2 The fractal dimension 

We have measured the fractal dimension in a number of independent ways^, and var- 
ious measurements agree. Here we limit ourself to report on one particular method, 
based on the distribution n N {r). From @ we have 

oo 

(r)i V = T7E'"Mr)~iV lAfe . (13) 

iV r=0 



Obviously, (13) could itself serve as a natural definition of du- By measuring n^j{r) 
we can record (r) N as a function of N and hence determine du- Strictly speaking, 
we expect this relation to be valid in the limit of infinite N, while finite N effects 
will be present in (|13|). The finite N effects should be dictated by the "the number 
of points" L corresponding to the linear size of the system, i.e. we expect 

(r) v const. const. . , 

' - const. + — - — + — - — + ■ • • . (14) 



ATi/dn L L 2 



2 The history of the seminal paper by Savvidy et.al. is interesting: The paper was rejected by 
4 computer journals, including Comput. Phys. Commun. where M. Liischer finally published a 
related paper and F. James the FORTRAN code of RANLUX in 1994. These ideas were commu- 
nicated to F. James by G.K. Savvidy during his stay at CERN in 1987. 

3 The details of these measurements will be reported elsewhere M. 
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Figure 1: The functions R a , N {d) for N = 2K, AK, 8K, 
a = 0.130 determined by minimizing Eq. ([T3) . 



, 25QK and for the optimal 



If we use the fact that N 1 ^ 11 is a typical measure for the linear extension of the 
manifold, Eq. ([14]) can be written, 

(r + ^^N^+Oi^), (15) 

by identifying L = N l l dH . The parameter a, which is considered the shift in r, incor- 
porates the next order correction. We will discuss possible physical interpretations 
of the shift a in detail in M. Now, let us define 

= ^^- < 16 > 

We determine the value of a and du in the following way: first we measure (r) N . 
for a certain number of different volumes N of the universes, ranging from N = 2K 
to N = 256K. For a given a we choose, for each couple N, Nj of N's, the dt^ such 
that 

Ra, Ni {d%) =Ra, Nj {4)- (17) 

For this choice of Ni, Nj we bin the data and estimate an error 5d%. Then we 
determine the average 

dH = ^—Y,d%, (18) 
# pairs ^ 



and compute 



(19) 
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2K-256K 
16K-256K 



3.57 



3.575 



3.58 



Figure 2: x 2 ( a ); defined by Eq. (|19"D for two sets of iVj's: Ni 
(solid curve) and Ni = 16K, 32K, . . . , 256K (dashed curve). 



2K, AK,..., 256K 



The preferred pair (a, G?ff (a)) is determined by the minimum of x 2 ( a )- This method 
works quite impressively. In Fig. [I] we have shown the intersection of the curves 
Ra,N{d) as a function of d for the optimal choice of a. The important point is that 
there exists a value of a where the curves intersect with high precision and that the 
the range of a where x 2 { a ) is acceptably small, i.e. 0(1), is quite small and hence 
dn will be determined with high precision. In Fig. |2| we show x 2 ( 
get 



0.139 ±0.005 



d H {a Ti 



3.574 ±0.003. 



In this way we 
(20) 



In (^) we have estimated the error as follows. Define an interval of acceptance 
[a min ,a max ] of a by demanding that x 2 ( a ) < %X* where \* — max{l, x 2 ( a m)} and 
find the variation of d(a) in this interval. After this we repeat the whole procedure 
by making various cuts in the pairs of AVs included in (|18D and (0), discarding 



successively the smallest AVs. The value of dn in Eq. (|20| ) agrees with the original 
value dii = 3.5 ± 0.2 ||. The strength of finite-size scaling is that one can obtain 
higher precision results with the use of much smaller lattices. The original simulation 
needed 50007^ size lattices. 



3.3 The boundary 

We now turn to the measurements of (l n ) r N . These observables are constructed from 
Pn(1,t), which can readily be measured in the simulations. If dim[iV] = dim[/ 2 ], 
then from scaling arguments, we expect 



r,N 



N n/2 F„ 



x 



X 



J\fl/d H 



(21) 
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Figure 3: The scaling functions F n (x) defined by Eq. (p2|) (left figure) and the scaling 
functions F n (x) defined by Eq. flSp, for n = 2 and N = 8K, 16K, . . . , 25QK. 



However, our measurements are consistent with the following scaling relation 

(l% N ~N 2n / d *F n (x), for n>2, (22) 
which implies that dim[Z] = dim[r 2 ]. Eq. (p2|) indicates that we have 

T)r,N ~ fOT 1 < T < iVVdH > " > 2 • ( 23 ) 

Again, relations like (^3) are expected to be valid up to finite size effects, as in 
Eq. flllj) . As a first phenomenological correction we use a shift r — > r + a as in 
(|I~5|) to find the best scaling function F n {x) for a suitable range of iVj's. In Fig. | 
we have shown F n (x) for n = 2 for the values of a which provide the best scaling 
function. Similar pictures exist for n = 3 and n = 4. This is to be compared to the 



scaling given in Eq. (21). We see that it is not possible to find a scaling function 
F n (x) if we used the ansatz (|2~T|). Assuming the scaling ( p2| ) we get an independent 
determination of dn = 3.63 ± 0.04. This result is remarkably consistent with the 
value of dn determined by the other methods we used, considering the systematic 
errors due to finite size effects. One, however, should consider the possibility of a 
dimensional relation of the form: 

dim[/] =dim[r 2 ( 1 ~ e )]. (24) 

In this case, for given dn, one can use the relation (l n ) rN = N 2n ( 1 ~ e ^ dH F n (x) in 
order to determine the value of e. Using the lowest value we obtained for d H using 
other methods, we obtain e < 0.03. Details will be published elsewhere 0. 

3.4 The spectral dimension 

Let us finally turn to the spectral dimension d s , as defined by fllPp. For a given 
triangulation /c7v(0,t) can be calculated by means of @ and (0), and the quantum 
average is then obtained by performing the average over different triangulations, 
since these are generated with the correct weight by the recursive sampling. In Fig. 
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Figure 4: A plot of In fcjv(0, t) versus lnt for N = 2K, 4K, . . . , 256K. 



|] we have shown a plot of In fcjv(0, t) versus Int. When t is too small, discretization 
effects interfere with the form (pUp, as discussed in details in ||. For larger values 



of t it is possible to perform a very good fit to (|10l), as is clear from Fig. [|. For the 
256K lattice the value of d s is 

d s = 1.980 ±0.014. (25) 

Due to finite size effects, the value of d s depends on the lattice size, but it clearly 
approaches 2 as N becomes large. Exploiting the N dependence of fcjv(r, t) it is 
possible to determine the exponents A and o introduced in Eqs. (PH|) and flnp. We 
find good agreement with the scaling prediction (|12"D . A detailed analysis of the 
diffusion equation and the numerical verification of (O) will be published elsewhere. 



4 Discussion 

The fractal structure of quantum gravity coupled to matter has provided us with 
somewhat of a puzzle. The first analytic result suggested that |T| 

d H = ~—, 7str = ^(c-l-v/(25-c)(l-c)). (26) 

Tstr J- * 

This formula has later been obtained in a number of different ways using the string- 
field approach developed in fl2| . There are several problems with the "proofs" . For 



instance, the string-field proof is based on the identification of proper time with 
geodesic distance. However, only in the case of pure two-dimensional gravity one 
can clearly identify the proper time used in string-field theory with the geodesic 
distance. Further, the formula ( ^6| ) predicts that d H — > oo for c — > 1 and d H ^ for 
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c — > — oo. The latter result is clearly undesirable since one expects that dn — ► 2 in 
the semi-classical limit c — ► — oo. In fact (^) predicts du = 2 in the case considered 
here, c = — 2 (7 str = —1), and a very recent constructive proof, using an explicitly 
constructed transfer matrix for c = —2, confirms this prediction [HJ. The prediction 
is clearly in disagreement with the computer measurements reported in this article. 
In a similar way the prediction (BQ) for c > has been in contradiction with the 
Monte Carlo simulations performed for the Ising model (c = 1/2, 7 str = —1/3) 
and the three-states Potts model (c = 4/5, 7 str = —1/5). In these cases one could 
argue that since the fractal dimension predicted is so large, the systems used in the 
computer simulations have been too small to observe the correct fractal dimensions. 
Although it is hard to understand how one can get all predicted i^Pi^exponents 
correct in the numerical simulations, but not being able to measure the fractal 
dimension, one could not rule out entirely this criticism. The important point is 
that the criticism is not valid for the present numerical investigation, since the 
predicted djj is small (namely 2), and we are able to deal with very large systems. 
The alternative prediction for dn is 

4 = -2^ = 2xf^ + f^. (27) 

The origin of this equation is to be found in the analysis of the diffusion equation 
in Liouville theory || and is based on the observation that (|8|)-([T2]) imply that 

dim[(r 2 (0) ] = dim[N 2 / dH ]. (28) 

This follows from the assumed scaling, independent of the specific model of dynam- 
ical triangulations. On the other hand, in Liouville theory one can use the De-Witt 
short distance expansion of the heat kernel in terms of geodesic distance to deduce 
that 

dim[(r 2 (t)) ] = dim[AH a - l/ai) ] , (29) 

provided r, N and t are viewed as continuum geodesic distance, continuum area and 
continuum diffusion time in Liouville theory. In ( 29[) a_ n denotes the gravitational 
dressing of a (n + 1, n + 1) conformal field, i.e. 



J d 2 ^^> n+l {g) - J d 2 ^g e a -^ n+1 (g) , for = e*® , 

where g^iO is the background metric. The requirement that e a -"^^ n+ i(g) is a 
(1,1) conformal field fixes 

In 

a n = . (30) 

1 + J(25-c-24n)/(25-c) 



For c = one obtains dn = 4 (in agreement with the transfer matrix prediction), 
while for c = —2 the prediction is 

d H ( c = -2) = (3 + Vl7)/2 = 3.561 . . . . (31) 
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It agrees with the numerical results reported in this paper. 

However, clearly this is not the complete story since we have also obtained 



dim[ l n } = dim[ r 2n ] , for n > 1 . (32) 

The same result is valid for c = and from numerical simulations for < c < 1 it 
seems to be valid in this region, too [18]. The reason such a result can appear is 
apparent from the transfer matrix calculations for c = 0. In this case we have 

p N (r, I) = ^G(l/r 2 ) + c 2 £- 3 /V(2 - l/r 2 )6(l - VZe) , (33) 

where e is a cut-off (the lattice spacing) and 

G{z) = (V 5 / 2 + -z~ 3 / 2 + —z 1 ' 2 ) e~ z , for N -> oo . (34) 

For low moments, n = and 1, the terms in ([J3|) which are singular for I —> will 
dominate the evaluation of 

(Hr,jv = ldU n p N (r,l), (35) 

while for n > 2 these terms become integrable and the cut-off dependence vanishes. 
It seems that the situation is the same at least for c G [—2, 1). The exists a regular 
part, (r,l), of pjv(r, Z) such that 

dlp% cg) (r, I) = d(l/r 2 )G(l/r 2 ) (36) 

is a function only of l/r 2 , and which dominates the integral (|35|) for n > 2, while a 
part, singular for I — > 0, dominates (p5|) for n = and 1. 

In order to fully understand the concept of fractal dimension we still have to 
provide an explanation of an expression like (133) for 0. 
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